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Eusocial Hymenoptera, such as the 
European honey bee, Apis mellifera, 
have the highest recombination rates 
of multicellular animals. 1 Recently, we 
showed 2 that a side-effect of recombina- 
tion in the honey bee, GC biased gene 
conversion (bGC), helps maintain the 
unusual bimodal GC-content distribu- 
tion of the bee genome by increasing 
GC-content in high recombination areas, 
while low recombination areas are losing 
GC-content because of biased AT muta- 
tions and low rates of bGC. Although 
the very high recombination rate of A, 
mellifera makes GC-content evolution 
easier to study, the pattern is consistent 
with results found in many other spe- 
cies including mammals and yeast. 3 Also 
consistent across phyla is the association 
of higher genetic diversity and divergence 
with high GC and high recombination 
areas. 4 '' Finally, we showed that genes 
overexpressed in the brains of workers 
cluster in GC-rich genomic areas with 
the highest rates of recombination and 
molecular evolution. 2 In this Addendum 
we present a conceptual model of how 
eusociality and high recombination rates 
may co-evolve. 

Recombination is high in several species 
with independent origins of eusocial- 
ity. 1 . 6 ' 7 j{ numDer 0 f studies have demon- 
strated fitness benefits in eusocial insects 
due to increased genetic variation in the 
colony, 8 " 10 and it was hypothesized that 
high recombination rates in eusocial 
insects were selected for because recombi- 
nation increases worker genotypic diver- 
sity. 1 ' 11 However, a recent modeling study 
demonstrated that recombination has 
little short-term impact on within-colony 
variation compared with factors such as 



polyandry. 12 What conditions favor the 
evolution of high recombination rates in 
social insects? Theoretical and empirical 
work points to two major factors: small 
effective population size, and strong selec- 
tion on linked genes. 13 ' 14 

Consider the early stages of the evo- 
lution of eusociality. First, the evolution 
of sociality is expected to greatly reduce 
effective population size 15 because of the 
larger amount of resources required to 
raise both workers and reproductives, 
higher variance in reproductive output 
between colonies and worker-produced 
males. 16 " 18 The reduction in effective pop- 
ulation size associated with the evolution 
of eusociality is expected to both increase 
linkage disequilibrium 19 and increase the 
frequency of slightly deleterious mutations 
via drift. 20 When recombination rate (R) 
is low, new advantageous mutations are 
less likely to fix if nearby deleterious alleles 
hitchhike along with them 21 — a phenom- 
enon called Hill-Robertson interference. 22 
Low R will also slow selection on advanta- 
geous alleles that are physically proximate 
but on different haplotypes, as the muta- 
tion with the highest selection coefficient 
will tend to drive the others to extinction 
before recombination can place them on 
the same haplotype. 23 The initial stages 
of the evolution of eusociality can thus 
be characterized by lower Ne, high link- 
age disequilibrium, higher genetic loads 
and reduced effectiveness of natural selec- 
tion due to interference between linked 
mutations. 

The above arguments suggest a para- 
doxical conclusion; the transition to euso- 
ciality created sub-optimal conditions 
for the evolution of caste-specialization 
and worker evolution via natural selec- 
tion. However, the concurrent evolution 
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of high recombination rates with euso- 
ciality — as observed across a number 
of independently derived eusocial lin- 
eages 1 — can reduce interference between 
linked mutations, thereby greatly enhanc- 
ing the effectiveness of natural selection. 
Assuming heritable variation in recom- 
bination rates, 24,25 queens with high rates 
of recombination are expected to produce 
more recombinant gametes allowing nat- 
ural selection to act more independently 
on physically linked mutations. This 
is expected to increase both queen and 
worker fitness by reducing interference 
and enhancing natural selection, which, 
in turn, should further select for higher 
rates of recombination. 

We propose three non-mutually exclu- 
sive processes that can exert strong selec- 
tive pressures for high recombination rates 
in eusocial insects. 

Divergence of Queen 
and Worker Phenotypes 

Caste specialization and division of labor 
is one of the hallmarks of eusociality. The 
genomes of social insects must generate 
two (or more) female phenotypes from 
the same genome (i.e., queens and work- 
ers). Several non-mutually exclusive pro- 
cesses can facilitate phenotypic divergence 
of castes, including caste-specific gene 
expression, alternative splicing, methyla- 
tion and gene duplication. 26 " 28 However, 
the end-result of the above processes is 
often very similar: mutations affecting 
worker phenotypes are likely to be physi- 
cally proximate to mutations affecting 
queen phenotypes. This will be true in 
cases where the same genes are expressed 
in both queens and workers but at differ- 
ent rates or at different times, 25 " 31 or in 
cases with caste-specific splice forms 32 and 
or caste-specific methylation patterns, 32 or 
if gene duplication events lead to physi- 
cally proximate caste-specific genes. In 
the above cases, we expect a general tug- 
of-war between linked mutations influ- 
encing worker and queen traits. 33,34 How 
is this tug-of-war resolved? Recent theory 
suggests that, all other effects being equal, 
queens invariably win because selection 
acts on them directly, but workers only 
experience indirect selection. 35,36 High 
recombination rates are expected to reduce 
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the tug-of-war between queens and work- 
ers, thereby allowing natural selection to 
optimize caste phenotypes (Fig. 1). 

Evolution of Worker Behavior 

Workers are the most abundant caste 
in eusocial insects, and their behavior 
strongly affects colony fitness. 37 By con- 
trasting the behavior and degree of caste 
specialization between primitively and 
advanced eusocial insects, it is clearly evi- 
dent how the evolution of worker behav- 
ior and specialization is associated with 
increased colony sizes and social complex- 
ity. 38 Given the major contributions of 
worker behavior to colony fitness, 8,39 it is 
inescapable that selection must be acting 
on genes that influence worker behavior 
and enhance colony fitness, 40 increasing 
both linkage disequilibrium and the odds 
of interference between linked mutations. 
Indeed, our recent work on a small subset 
of honey bee genes associated with worker 
behavior has shown a prevalence of natu- 
ral selection and high linkage disequilib- 
rium near selective sweeps. 41 We have also 
demonstrated that genes associated with 
behavior and neurobiology were ances- 
trally found in regions of the genome with 
above average GC-content. 2 GC-rich areas 
are expected to contain more mutations 2 
and recombination motifs. 42 Adaptive evo- 
lution of worker behavior is therefore likely 
to increase linkage disequilibrium and 
interference selection in GC-rich regions, 
which in turn can select for the evolution 
of additional recombination hotspots. 42 

Indirect Genetic Effects 
between Workers 

Although selection on mutations affect- 
ing worker traits is relatively weak because 
it is indirect (i.e., as it changes the fitness 
of their reproductive kin), 35 social inter- 
actions between workers can cause mea- 
surable indirect effects on colony fitness 43 
leading to social regulatory networks 
shaped by colony-level selection. 44 Indirect 
phenotypic effects occur when individuals 
change phenotypes as a result of interac- 
tions with others. 45 Recent theory 46 sug- 
gests that behavioral interaction effects 
may strongly increase the selection coeffi- 
cients of mutations influencing interaction 
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traits; this increase in selection coefficients 
is in theory proportional to the number of 
interactions among individuals 46 and is 
expected to be greater in larger colonies. 
If colony fitness increases due to increased 
worker-worker interactions, selection on 
mutations with indirect effects on workers 
can be substantial. Similar to scenario 2, 
selection on mutations with indirect 
genetic effects on worker behavior will 
increase linkage disequilibrium and inter- 
ference selection, which in turn can select 
for the evolution of higher recombination 
rates. 

We propose that the evolution of euso- 
ciality and recombination rates have been 
linked in a feedback loop that has strongly 
affected genome evolution (Fig. 2) . Higher 
recombination rates were selected for dur- 
ing the early stages of the evolution of 
sociality because they reduce interference 
selection allowing natural selection to 
optimize caste phenotypes. This enhanced 
division of labor then allowed colony sizes 
to increase, which in turn drives additional 
reductions in effective population size and 
concurrent selection for higher recombina- 
tion rates. 

Some components of the loop cause 
neutral or non-adaptive patterns of molec- 
ular and genome evolution, such as the 
effects of reduced effective population size 
as colonies grow larger, and the effects of 
biased gene conversion on GC content. A 
side effect of increased GC-content due to 
biased gene conversion may be to increase 
the number of potential binding sites for 
fast-evolving recombination hotspot bind- 
ing factors (such as PRDM9 in mam- 
mals 47,48 ). This phenomenon, in turn, 
generates new recombination hotspots in 
GC-rich regions, which replaces recombi- 
nation hotspots that burn out. 48 Our model 
suggests that not all unusual genomic fea- 
tures found in eusocial insects were caused 
by adaptive evolutionary forces. 

Our model leads to several testable 
predictions. Our assumption that muta- 
tions affecting queens and worker traits 
are physically proximate — representing 
one of three possible scenarios that can 
select for high R — can be tested with 
more knowledge on the pattern of gene 
expression, methylation, alternative splic- 
ing and gene duplication in sequenced 
eusocial genomes. Emerging techniques, 
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such as genome wide association mapping, 
will allow us to link genetic polymor- 
phisms with queen and worker traits in 
the foreseeable future. Further, population 
genomic studies will inform us about the 
prevalence of natural selection and its tar- 
gets in eusocial organisms. Additionally, 
our model predicts a negative relationship 
between recombination rates and GC con- 
tent with effective population size across 
social species. Finally, our model predicts 
some fundamental differences in the 
degree of linkage disequilibrium, interfer- 
ence selection and the efficiency of positive 
selection, between the early and late stages 
of the evolution of sociality. It would be 
interesting to contrast advanced eusocial 
and primitively eusocial taxa, and their 
closely related solitary ancestors, at some 
of the key population genetic parameters 
discussed herein. Studies on species that 
recently reverted from social to solitary 
lifestyles 49 are also relevant, as we believe 
that these species may not have passed the 
evolutionary hurdle imposed by interfer- 
ence selection during the initial stages of 
the evolution of sociality. 

We hope our conceptual model stimu- 
lates work on genomic and population 
genomic studies of eusocial taxa, as well as 
theoretical studies on the extent to which 
haplodiploidy and eusociality provide the 
necessary conditions for the evolution of 
high recombination rates. 
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Figure 1. Recombination and the evolution of caste specialization. Recombination can enhance 
concurrent selection on initially linked mutations that influence queen and worker traits. We 
assume that both mutations are recent (i.e., have low frequency). We depict frequencies of the 
different gametes assuming a recombination rate r, which ranges between 0 to 50% (i.e., free 
recombination). (A) A haplotype containing a slightly deleterious mutation affecting queens 
linked to a beneficial mutation affecting worker traits and enhancing colony fitness. Recombi- 
nation can generate a haplotype containing the beneficial worker mutation and the wild-type 
queen allele, thereby facilitating both positive selection (on the 'worker' mutation) and purifying 
selection (on the "queen" mutation). (B) A haplotype containing a beneficial mutation affecting a 
queen trait linked to a slightly deleterious mutation affecting a worker trait. Recombination can 
facilitate concurrent positive selection (on the 'queen' mutation) and purifying selection (on the 
"worker" mutation). (C) Two haplotypes, each containing a beneficial mutation affecting a queen, 
and a worker trait respectively. These two mutations are on different haplotypes but are spatially 
proximate. Without recombination, selection will fix one haplotype at the expense of the other. 
However, one of the recombinant gametes will contain both beneficial mutations allowing for 
concurrent positive selection on queen and worker traits. 
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Figure 2. Conceptual model for the evolution of eusociality and recombination rates. Increasing 
social organization and caste complexity lead to colonies with more workers per queen, reduc- 
ing effective population size N e . Lower N e increases linkage disequilibrium and fixation of mildly 
deleterious mutations, especially in genomic regions with low recombination. This increases Hill- 
Robertson interference that can in turn select for higher R, which allows natural selection to act 
more efficiently on mutations affecting worker traits such as behavior. This in turn facilitates the 
evolution of higher social complexity and colony sizes that further depress N g . 
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